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Abstract 

I discuss methods of calculation of propagator diagrams (massless, those of Heavy 
Quark Effective Theory, and massive on-shell diagrams) up to 3 loops. Integration- 
by-parts recurrence relations are used to reduce them to linear combinations of basis 
integrals. Non-trivial basis integrals have to be calculated by some other method, 
e.g., using Gegenbauer polynomial technique. Many of them are expressed via hyper- 
geometric functions; in the massless and HQET cases, their indices tend to integers 
at e — >• 0. I discuss the algorithm of their expansion in e, in terms of multiple C 
values. These lectures were given at Calc-03 school, Dubna, 14-20 June 2003. 
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1 Introduction 



We shall discuss methods of calculation of propagator diagrams in various theories and 
kinematics, up to 3 loops. In order to obtain renormalization-group functions, it is often 
enough to calculate propagator diagrams. Massless diagrams are discussed in Sect. |21 
they are useful, e.g., in QCD with light quarks. Heavy Quark Effective Theory (HQET) is 
discussed in Sect. El and massive on-shell diagrams - in Sect.lU Finally, some mathematical 
methods, most useful in massless and HQET calculations, are considered in Sect. El 

Throughout these lectures, we shall discuss scalar Feynman integrals. Tensor integrals 
can be expanded in tensor structures, scalar coefficients are found by solving linear sys- 
tems. Their solution can be written as projectors (tensor, 7-matrix) applied to the original 
diagram. 

Scalar products of momenta in the numerator can be expressed via the denominators. 
In some classes of problems, there are not enough independent denominators to express 
all scalar products. Then we add the required number of scalar products to the basis, and 
consider scalar integrals with powers of these scalar products in the numerators. 

Some diagrams contain insertions into internal lines. The denominators of the prop- 
agators on both sides coincide. If powers of these denominators are rii and 77-2, we can 
combine them into a single line with the power rii + n2 (Fig. QJ. 



Throughout these lectures, Minkowski notations with the metric + — — — are used. 




ni + 77-2 



Figure 1: Insertion into a propagator 



2 Massless propagator diagrams 



2.1 1 loop 

The 1-loop massless propagator diagram (Fig. ^ is 




(2.1) 



-iO. 



It is symmetric with respect to 1 ^ 2. 
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k + p 



Figure 2: 1-loop massless propagator diagram 

If rii is integer and ni < 0, it becomes the vacuum massless diagram (Fig. E> with the 
polynomial in the numerator. It reduces to a sum of powers of — multiplied by 

integrals 



/ 



0, (2,2) 



This integral has the mass dimension d — 2n; it contains no dimensionful parameters, and 
the only possible result for it is 0^. Therefore, G{ni,n2) = for non-positive integer rii (or 
71,2). More generally, any massless vacuum diagram with any number of loops is scale- free 
and hence vanishes. 




Figure 3: 1-loop massless vacuum diagram 
This diagram is calculated by the Fourier transform to x-space and back: 



;-p2 - iO)n (^27rY T{n) {-x^ + ioy/^-^ 

Our diagram in x-space is the product of two propagators ()2.3p with the powers rti and 77-2: 



_^T{d/2-ni)T{d/2-n2) 1 

T{ni)T{n2) (_^2-)d-ni-n2 



^This argument fails at n ~ (i/2; more careful investigation ^ shows that the right-hand side contains 
5{d~2n). 
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Transforming it back to p-space ()2.4|) . we arrive at 

r{-d/2 + ni + n2)T{d/2 - ni)r{d/2 - ns) 



G{ni,n2) 



T{ni)T{n2)T{d - rii - 722) 



Diagrams are calculated by going to Euclidean momentum space: ko = ikEo, k"^ 



(2.5) 



Using a dimensionless loop momentum K = ks/ a/— p^, we can rewrite ()2.H) as 
d'^K 



where n is a unit Euclidean vector (n^ = !)• We can perform inversion 



K 



K' 



K'2 



d^K 



d'^K' 



Then 



{K + n 



and we obtain (Fig. E} 

G(ni, 722) = G(ni, d — rii — 722) . 



^2 




d — rii — n2 




rii 



(2.6) 



(2.7) 



Figure 4: Inversion relation 



2.2 2 loops 

There is 1 generic topology of 2- loop massless propagator diagrams (Fig. E}: 

"^'^^ "^'^^ -A-P'r^ "^G(ni, ^2, ^3, n4, n,) , (2.8) 



Di = -{ki+p)'\ D2 = -{k2+p)\ Ds = -kl, D^ = -kl, D5 = -ik,-k'^ 
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All other diagrams (e.g., Fig. are particular cases of this one, when some line shrinks to 
a point, i.e., its index rij = 0. This diagram is symmetric with respect to (1 ^ 2,3 ^ 4), 
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Figure 5: 2- loop massless propagator diagram 



and with respect to (1 ^ 3,2 ^ 4). If indices of any two adjacent lines are non-positive, 
the diagram contains a scale-free vacuum subdiagram, and hence vanishes. 
If ns = 0, our diagram is just the product of two 1-loop diagrams (Fig. in}: 



G{ni, 712, ^3, n4, 0) = G{ni,n3)G{n2, n^) . 
If Til = 0, the inner loop gives G{ns, n^) / {—k'^)"'^'^^^''^^'^ , and (Fig. Hj) 

G(0, n2, ris, n^, n^) = G^ris, n^)G{n2, Ui + n^ + n^- d/2) . 
The cases n2 = 0, ns = 0, 71,4 = are symmetric. 




(2.9) 



(2.10) 



n2 




Figure 6: Trivial case ns = 




n^ + ris + nc,- d/2 




715 



n2 



n2 



Figure 7: Trivial case ni = 

When all > 0, the problem does not immediately reduce to a repeated use of the 
1-loop formula ()2.5|) . We shall use a powerful method called integration by parts j2]. 
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It is based on the simple observation that any integral of d/dki{- ■ ■) (or d/dk2{- ■ ■)) 
vanishes (in dimensional regularization no surface terms can appear). From this, we can 
obtain recurrence relations which involve (^(rai, n2, ns, 71,4, ns) with different sets of indices. 
Applying these relations in a carefully chosen order, we can reduce any G{ni,n2, n^, n^, n^) 
to trivial ones, like ^TT^ . 

The differential operator d/dk2 applied to the integrand of ()2.8|) acts as 

^ ^2{k2 + p) + ^2k2 + ^2{k2 - k,) . (2.11) 



dk2 D2 ' ' A D5 

Applying {d/dk2) ■ k2 to the integrand of ()2.8|) . we get a vanishing integral. On the other 
hand, from (|2SD, 2k2 ■ k2 = -2^4, 2{k2 + p) ■ ^2 = (V) - ^2 - A, 2(^2 - k^) ■ k2 = 
D3 — D4 — D5, we see that this differential operator is equivalent to inserting 

d-n2-n,-2n, + ^(( V) - A) + ^(^3 - A) 
^2 ^5 

under the integral sign (here {d/dk2) ■ k2 = d). Taking into account the definition ()2.8j) . 
we obtain the recurrence relation 

- ns - ns - 2n4 + n22+(l - 4") + n55+(3~ - 4")] G = . (2.12) 

Here 

l^G(ni, n2, ^3, n4, n^) = G{ni ± 1, 77,2, ng, n4, 77,5) , (2.13) 

and similar ones for the other indices. 

This is a particular example of the triangle relation. We differentiate in the loop 
momentum running along the triangle 254, and insert the momentum of the line 4 in the 
numerator. The differentiation raises the degree of one of the denominators 2, 5, 4. In the 
case of the line 4, we get —2D4 in the numerator, giving just — 2n4. In the case of the line 
5, we get the denominator of the line attached to the vertex 45 of our triangle, minus 
the denominators D4 and D^. The case of the line 2 is similar; the denominator of the line 
attached to the vertex 24 of our triangle is just — p^, and it does not influence any index 



of G. Of course, there are three more relations obtained from ()2.12|1 by the symmetry. 
Another useful triangle relation is derived by applying the operator {d/dk2) ■ {k2 — ki): 

[d-n2-n^- 2^5 + n22+(l" - 5") + n44+(3" - 5")] G = 0. (2.14) 

One more is obtained by the symmetry. Relations of this kind can be written for any 
diagram having a triangle in it, when at least two vertices of the triangle each have only a 
single line (not belonging to the triangle) attached. 

We can obtain a relation from homogeneity of the integral ()2.8j) in p. Applying the 
operator p- (d/dp) to the integral ()2.8|) . we see that it is equivalent to the factor 2{d — J2^i)- 
On the other hand, explicit differentiation of the integrand gives —{ni/Di){—p'^ + Di — 
Ds) - {n2/D2){-p^ + D2- D4). Therefore, 

[2{d - ns - 714 - ^5) - ni - 712 + nil+(l - 3") + 7222+(l - 4")] 1 = 0. (2.15) 
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This is nothing but the sum of the {d/dk2) ■ /c2 relation (j2.12j) and its mirror-symmetric 
[d/dkij ■ ki relation. 

Another interesting relation is obtained by inserting {ki+p)^ into the integrand of ()2.8|l 
and taking derivative d/dp^ of the integral. On the one hand, the vector integral must be 
proportional to p^, and we can make the substitution 

{ki+p)-p f Di- p 
ki+p^ 5 p - ' ' ' 



in the integrand. Taking d/dp^ of this vector integral produces ()2.8|) with 



^ ^-^ / \ — 

inserted into the integrand. On the other hand, explicit differentiation in p gives 

d + ^2ik^ + pf + ^2ik2 + p) ■ {k, + p) , 

JJl U2 

2ik2+p)-ik,+p) = D,-D,-D2. 
Therefore, we obtain 

ici + ni-n3-n4-n5+(|rf-^ni)(l~-3") + n22+(l"-5-) G = 0. (2.16) 

This relation has been derived by S.A. Larin in his M. Sc. thesis. Three more relations 
follow from the symmetries. 

Expressing G(ni, 77-2, "n-s, n4, 77,5) from ()2.14|) : 

N n22+(5- - 1-) +n44+(5- - 3-)^ 
G ni,n2,n3,n4,n5 = ^ -G , 2.17 

we see that the sum rii + ns + reduces by 1. If we start from an integral belonging 
to the plane ni + + = n in Fig. |Hl then each of the integrals in the right-hand side 
belong to the plane rii + + = n — 1 (Fig. IS)). Therefore, applying ()2.17|) sufficiently 
many times, we can reduce an arbitrary G integral with integer indices to a combination 
of integrals with 71,5 = (Fig. IHl ()2.9|) ). rzi = (Fig. [71 ()2.10p ). 723 = (mirror-symmetric 
to the previous case). Of course, if max(7T,2, 7^4) < max(?7,i, 723), it may be more efficient to 
use the relation mirror-symmetric to ()2.14j) . The relation ()2.1(i|l also can be used instead 
of (EH- 

Let's summarize. There is 1 generic topology of 2-loop massless propagator diagrams 
(Fig.ini). The integrals with all integer rti can be expressed as linear combinations of 2 basis 
integrals (Fig. ITIHl : coefficients are rational functions of d. Here G„ are 77,- loop massless 
sunset diagrams (Fig. HH): 

n ^ I T{l + ne)T-+\l-e) 

" (n + l-7zf)J(n + l)f-2n-l)^ V{l-{n + l)e) ' ' 
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Figure 8: Reduction oi n = ni + + 




Figure 9: 2- loop massless propagator diagram 




Figure 10: Basis diagrams (all indices equal 1) 




Figure 11: n-loop massless sunset diagrams 
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Using inversion ()2.6|) of both loop momenta and 



l\2 



we obtain relation of Fig. IT^ . 




d — rii — n 



— n 




Figure 12: Inversion relation 



2.3 3 loops 

There are 3 generic topologies of 3-loop massless propagator diagrams (Fig. IT!?|l . Each 
has 8 denominators. There are 9 scalar products of 3 loop momenta ki and the exter- 
nal momentum p. Therefore, for each topology, all scalar products in the numerator can 
be expressed via the denominators and one selected scalar product. Integration- by-parts 
recurrence relations for these diagrams have been investigated in [2]. They can be used 
to reduce all integrals of Fig. El with arbitrary integer powers of denominators and ar- 
bitrary (non-negative) powers of the selected scalar product in the numerators, to linear 
combinations of 6 basis integrals (Fig. IT^ . This algorithm has been implemented in the 
SCHOONSCHIP package MINCER [3 and later re-implemented [H in FORM It has 
also been implemented in the REDUCE ^01 package Slicer ^2] Only the last, non-planar, 
topology in Fig. IT3l involve the last, non-planar, basis integral in Fig. El 

The first 4 basis integrals are trivial: they are expressed via Gn ()2.18j) (Fig. ITT|l . and 
hence via F functions. The 4-th one differs from the 3-rd one (G3) by replacing the 2- 
loop subdiagram: the second one in Fig. HUl (G'2/(— A;^)'^"'^) by the first one (G'^/(— A;^)^"'^). 
Therefore, it can be obtained from by multiplying by 

GlG{lA-d) _2d-hG\ 
^2^(1, 3 -rf) ~ d-'iOj' 

^This is one of the elements of the symmetry group ^2 x S'e El (with 1440 elements) of the diagram of 
Fig-El If we connect the external vertices by a line having hq ~ X]i=i obtain a logarithmically 

divergent 3-loop tetrahedron vacuum diagram It is proportional to our original diagram times the logarithm 
of the ultraviolet cut-off. We can cut it at any of its lines, and get a diagram of Fig. El with 5 indices 7ii out 
of 6. This gives the tetrahedron symmetry group. Also, we have the inversion relations fFig. I12|. Fourier 
transform to x-space gives a new integral of the form of Fig. with rii — > | — n^. There is an additional 
symmetry following from the star-triangle relation ^jSl. 
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Figure 13: Topologies of 3- loop massless propagator diagrams 




rig d — uq — rij — ris 




Figure 15: Inversion relations 




Figure 16: Inversion relation: all Ui = 1, d = A 
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The 5-th basis integral is proportional to the 2-loop diagram G(l, 1, 1, l,n) with a non- 
integer index of the middle line n = e, and will be discussed in Sect. 12.41 The 6-th one, 
non-planar, is truly 3-loop and most difficult; it will be discussed in Sect. 12.51 

Performing inversion ()2.(ij) of the loop momenta, we obtain the relations in Fig. ^1 For 
example, the ladder diagram with all indices = 1 is convergent; its value at = 4 
is related to a simpler diagram (Fig. [TB|) by the second inversion relation. The non- 
planar topology (Fig. [T^ involves lines with sums of 3 momenta; they don't transform 
into anything reasonable under inversion, and there is no inversion relation for it. 



2.4 G(l,l,l,l,n) 
This diagram is ^3] 

G'(l,l,l,l,n) = 2r(f-l)r(f 

2r (i - 1) 

{d-2n- A)r{n + 1)1 (|d - n - 4) 



n 



- 1) T{n 

3-^2 



d + 3) X 



1, (i — 2, n 
n + 1, n - 



+ 2 



(2.19) 

TT cot 7i{d — n) 
T{d-2) 



This is a particular case of a more general result ^Hj for the diagram with 3 non-integer 
indices. The two lines with unit indices must be adjacent (Fig. [T2j); all such diagrams are 
equivalint, due to the tetrahedron symmetry mentioned above. This result was derived 
using Gegenbauer polynomial technique ^1]; we shall discuss it for a simpler example in 
Sect. 13.61 When n is integer, the 3F2 in ()2.19|) is expressed via F functions; in order to 
show equivalence to the standard results (Sect. 12. 2j) . one should use 



T{l + le)T{l-le) 



Tile 



sin -nle 

and trigonometric identities. 
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^2 




^3 



n2 



Figure 17: Diagrams calculated in jT3] 

Similar expressions for these diagrams were also obtained in ^H]^- In our particular 

■^The statement that these diagrams can be expressed via 3_F2 functions of unit argument was commu- 
nicated by D. Broadhurst to A. Kotikov before |13|. 
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case, they reduce to 

{d -3){d- 4)r(n)r (frf - n - 4) 



2T{n + 3-d)T^ (|_l)r(f 



n 



1) 



^(1,1,1,1,^) 



l,j-2,n + 3-d 



M~2n- 10 
d-n-S^^V n,n + A-d 

+ / T^TTCotTT n - - 2r I - 1) 

r(c/-4)r(f -1) ^ ^ \2 ) 

3d-2n- 10 
d — n — 3 



"3-^2 



lA — n,d — n — 3 



2' 



d — n 



1 



+ 



r(n)r (|c/-n-4) 

r(rf-4)r(f-i) 
1 



TT cot TT 



d-4 



n 



- 2n - 2 



3-^2 



^ d — TT, — 3 

3j „ I 1 d 
2 



1, n + 5 — |d, n + 1 2 
3-f,n + 2-f 



r(f-i) 



(2.20) 



+ 



-4 



T{n)T {\d-n-A) 

r(d-4)r(f -1) 

X /id 



;rcot7r^-2— ^-A_r(|-1) 



d-2n-2 



3-^2 



'2 "'2 

frf - n - 4, 1 - n 



d - 2n - 2 
1 



+ 



T{n)T (|rf-n-4) 

r(d-4)r(f-i) 



TT 



cot7r(f-n)-2r(f-l) . 



Kazakov jH] obtained an expression for G{1, 1, 1, l,n) via two 3F2 functions of argument 
— 1 much earher. Recently, G(ni, n2, na, 77,4, 77,5) for arbitrary indices was calculated ^Hl in 
terms of a double integral or several double series. These series can be systematically ex- 
panded in e up to any order, and coefficients are expressed via multiple C- values fSect. l^TTj) . 

In order to calculate the 2- loop diagram of Fig. IHl with a 1-loop insertion in the middle 
light line, we need G(l, 1, l,l,n + e). It is easy to shift by ±1 using the relation [2j 

[{d - 2^5 - 4)5+ + 2{d - ris - 3)] ^(1, 1, 1, 1, n^) = 2 l+(3- - 2-5+)G'(l, 1, 1, 1, n,) , 

(2.21) 

which follows from integration-by-parts relations (Sect. I2.2|l (all terms in its right-hand 
side are trivial). Therefore, it is sufficient to find it just for one n. The simplest choice for 
which the algorithm of e-expansion (Sect. Ej) works for (|2.19|) is ns = 2 -|- e: 



, 2r(l + 3e)r(l -e) 
G(l, 1, 1, 1, 2 + e) = ^ X 



1+25 



r(l - 2e)T{l-e) 



(2 + £)(l + £)2r(l-4£)r(l + £ 



"3-^2 



l,2-2£,2 + 2e 
3 + e,3 + 2e 



TT cot 3ne 
26{l-2e) 



(2.22) 
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1 =4(C2-l) + 6(-4C3 + 3C2-l)^ 



Expansion of this 3-F2 in e is (Sect. 

/ l,2-2e,2 + 2e 
^ ^\ 3 + e,3 + 2e 

+ 2(41(4 - 54C3 + 22C2 - 9)e^ + 3(-124C5 + 24(2(3 + 123(4 - 88(3 

Similarly, we obtain from ()2.2()j) for = 1 + e 

3e%l - 2e)T{l + e)T{l - 4e) 



30(2 - 8)e 



r2(i-e)r(i-2e)r(i + 3e) 



4 ^ 



1, —5, 35 
1 + e, 1 + 3e 



_1 

3 ^ ^ V 1 + e, 1 - 3e 



1, —e, —3e 



2 



l,45,2e 
1 + e, 1 + 2£ 



■G(l,l,l,l,l + £) 

r(i + £)r(i-4£) 
r(i -£)r(i -2e) 
r(i + £)r(i-4£) 



2 ^ H 1 - 4e, 1 - 2£ 



r(i 


-e)T{l 


-25) 


r(i 


+ s)V{l 


-As) 


r(i 


-em 


-2e) 


r(i- 


fe)r(l- 




r(i- 


-.)r(i- 


-2e) 



Tie cot 37re 



Tre: cot vre 



Tre: cot vre 

2 



where 

3-^2 

3-^2 
3-^2 
3-^2 



1, — £, 3£: 
1 + e, 1 + 3£: 

1, — e, — 35 
1 + e, 1 - 3£ 

l,45,2e 
l + e,l + 2e 

1,-e, -2e 
1 - 4e, 1 - 2e 



1 - 3(2^' + 18(36^ 
69 



123 



Tie cot 2716 



(46^ + 9(31(5 -6(2(3)5' + 



1 + 3C2e' 



^,^4 ^27(-3(5 + 2(2(3)5' + 



1 + 8(26^ + 92(46^ + 72(-3(5 + 2(2(3)^' + ■ • • 
= 1 + 2(2^' + 18(36^ + 101(46^ + 18(23(5 + 2(2(3)5' + 



Also, the last formula in (j2.2()|l gives for = e 
G(l, 1,1.1,.)- 2r(l-.)r(l+3e) 



3e(l - 2£)(1 - 3e)(l - 4e)r(l + e)T{l - Ae) 



X 



1 - 2e 



r(i-£)r(i- 25)3^2 



1,-e, 1 -2£ 
2 - 4£, 2 - 2£ 



+ (1 - 4e)r(l + e)T{l - Ae)Ti cot 27re - r(l - e)T{l - 2e) 



where 



3-^2 



l,-£, 1 -2£ 
2 - 4e, 2 - 2£ 



101 



+ (— (4-45(3 + 3(2-6)5^ 



1 + ((2 - 2)e + (9(3 - 5(2 - 3)e^ 



505 



(2.23) 

3 , 



(2.24) 



(2.25) 



(2.26) 



(2.27) 



207(5 + 18(2(3 - —(4 + 27(3 + 3(2 - 15 1 + 
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It is not difficult to find several additional terms to ()2.23p . ()2.25p . ()2.27p . Of course, all 
4 series (IT^ for G(l, 1, 1, 1, 1 + e) agree, and those for 1, 1, 1, 2 + e) ^FT^ and 
1, 1, l,e) ()2.26|) agree with the recurrence relation ()2.21|1 . 

2.5 Non-planar basis integral 

The non-planar basis diagram (the last one in Fig. IT^ is finite; its value at e = can be 
obtained without calculations by the method called gluing j2]. 

Let's consider the 4-loop vacuum diagram in the middle of Fig. ^1 Let all lines have 
mass m, then there are no infrared divergences. There are no divergent subdiagrams; the 
diagram has an overall ultraviolet divergence, and hence a 1/e pole. 




Figure 18: Gluing 

We can imagine that the middle vertex consists of 2 vertices connected by a shrunk line 
(with power 0), one vertex involves 2 lines on the left, and the other one - 2 lines on the 
right. If we cut this shrunk line, we get the ladder diagram (the upper one in Fig. ITH|l . It 
is finite. Therefore, the ultraviolet divergence of the vacuum diagram comes from the last 
integration in d'^p of the ladder diagram. At p ^ cxd, we may neglect m; by dimensionality, 
the massless ladder diagram behaves as {A + 0{e)/{—p'^)'^~^^'^ . Reducing it to the basis 
integrals and using results of Sect. 12. 4[ we can obtain A = 20^5 (first line of Fig. IT8|l "^. It is 
not difficult to obtain several additional terms. The ultraviolet divergence of the vacuum 

^this 0(1) term was found in j2j using Gegenbaucr polynomials in x-space |14) : however, it is difficult 
to obtain further terms of e expansion by this method. 
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diagram comes from 



df-p 



J (•_p2^2+3e 



UV 



9 X-Se 1 



r(2 -e) J r(2 - e) Ae 



where A is an infrared cutoff, see the middle hne in Fig. ^1 

On the other hand, we can imagine that one of the vertices in the middle of the vacuum 
diagram involves the upper left line and the lower right one, and the other vertex - the 
lower left line and the upper right one. Cutting the shrunk line, we get the non-planar 
diagram (the lower one in Fig. ITHj) . It is finite; at p — > oo, where m can be neglected, it 
behaves as {B + 0{e) / (— p^)^"*"^^. Integrating it in d'^p, we obtain the same vacuum diagram, 
with the same 1/e pole. Therefore, B = A: the non-planar basis diagram has the same 
value at d = A as the ladder one |2j, namely, 20^5 (Fig. I18|l. This method tells us nothing 
about further terms of expansion in e. The highly non-trivial problem of calculating the 
0{e) term in the non-planar basis diagram was solved in 



3 HQET propagator diagrams 
3.1 Crash course of HQET 

Heavy Quark Effective Theory (HQET, see, e.g., [IZl CHI 1131201) is an effective field theory 
constructed to approximate results of QCD for certain problems with a single heavy quark. 
Let the heavy quark stay approximately at rest (in some frame); its momentum and energy 
are \p\ < -E, \po — m\ ^ E. Momenta and energies of light quarks and gluons are \ki\ < -E, 
l^oil ^ E. Here E <^ m is some fixed characteristic momentum, and we consider the 
limit m oo. Scattering amplitudes and on-shell matrix elements of operators in QCD, 
expanded up to some order in E/m, can be reproduced from a simpler theory - HQET. 
The lowest-energy state of this theory ("vacuum") is a single heavy quark at rest, and 
has energy m. Therefore, for any system containing this quark it is natural to measure 
energy relative this zero level, i.e., to consider its residual energy po = Po — The 
free heavy quark has dispersion law (dependence of energy on momentum, or mass shell) 
Po = \/ 'w? + p^ — m. Neglecting 1/m corrections, we may simplify it to Po{p) = 0: the 
heavy-quark energy is zero, and does not depend on its momentum. The heavy quark at 
rest is described by a 2-component spinor field Q, or a 4-component spinor having only 
upper components: 70Q = Q- The dispersion law po{p) = follows from the Lagrangian 
L = Q^idoQ. Reintroducing the interaction with the gluon field by the requirement of 
gauge invariance, we obtain the leading-order (in 1/m) HQET Lagrangian 

L = Q+iD^Q + --- (3.1) 

where all light-field parts (denoted by dots) are exactly the same as in QCD. The heavy 
quark interacts with Aq: creates the coulomb chromoelectric field, and reacts to external 
chromoelectric fields. 
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The Lagrangian ()3.H) gives the heavy quark propagator 

5(p) = :^^, Six) = -teixo)6ix,). (3.2) 
Po + ^0 

In the momentum space it depends only on pq but not on p, because we have neglected the 
kinetic energy. Therefore, in the coordinate space the heavy quark does not move. The 
vertex is igv'^t"', where = (1,0,0,0) is the heavy-quark 4-velocity. Heavy-quark loops 
vanish, because it propagates only forward in time. 

HQET is not Lorentz- invariant, because it has a preferred frame - the heavy-quark 
rest frame. However, it can be rewritten in covariant notations. Momentum of any system 
containing the heavy quark is decomposed as 

p = mv + p , (3.3) 

where the residual momentum is small: \p^\ < E. The heavy-quark field obeys -^Q^ = Qy. 
The Lagrangian is 

L = Qyiv ■ DQy + • • • (3.4) 
It gives the propagator 
1 + / 1 

Sip) = ^—JL^ (3.5) 

and the vertex igv^t"-. 

The QCD propagator at large m becomes 

s(p) = 4±i!L = ^^(Lti)+f = i±^J_ + ofiV (3,6) 

_ jy^z 2mp - f + p^ / p ■ V \m 

A vertex ig'jH^ sandwiched between two projectors may be replaced by igv^t"". There- 
fore, at the tree level all QCD diagrams become the corresponding HQET diagrams, up 
to 1/m corrections. In loop diagrams, momenta can be arbitrarily large, and this corre- 
spondence breaks down. Renormalization properties of HQET differ from those of QCD. 
The ultraviolet behavior of an HQET diagram is determined by the region of loop mo- 
menta much larger than the characteristic momentum scale of the process but much 
less than the heavy quark mass m (which tends to infinity from the very beginning). It 
has nothing to do with the ultraviolet behavior of the corresponding QCD diagram, which 
is determined by the region of loop momenta much larger than m. 

3.2 1 loop 

The 1-loop HQET propagator diagram (Fig. ITT?|l is 



kp + UJ 2 

Di = , D2 = -k 



(3.7) 
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ko + uj 



Figure 19: 1-loop HQET propagator diagram 

It vanishes if tt-i < or 77-2 < 0. 

The Fourier transform of the 1-dimensional (HQET) propagator is 

' "^"^ ' -r-'e{t), (3. 



-uj - iO)" 27r r(n) 



-00 
00 



e t dt- ^_^_^Q)„+i • (^-JJ 







Our diagram in x-space is the product of the heavy propagator ()3.8|1 with the power ni 
and the hght propagator ()2.3|) with the power 77,2- The heavy quark stays at rest: x = 0. 
Therefore, —x"^ = [itY in ()2.3|) ^. Our diagram in the x-space is 

r(ni)r(n2) ' ^' 

Transforming it back to p-space ()3.9|1 . we arrive at 
jt ^ r(-rf + r2i + 2n2)r((i/2-n2) 

/ ni, ns = . 3.10 

r(r2i)r(n2) 

3.3 2 loops 

There are 2 generic topologies of 2-loop HQET propagator diagrams. The first one is 
(Fig. mi): 

-n\-2uf'^'~-^-^--^h{n,, n2, ns, n,, n,) , (3.11) 



Di = , D2 = , D-i = -k^ , Di = -k^ , D5 = -{ki - k2) . 

UJ UJ 

This diagram is symmetric with respect to (1 ^ 2,3 ^ 4). It vanishes if indices of any 
two adjacent lines are non-positive. 

^Why not {—ity^l The Wick rotation to the Eudidean x-space \s t = —He- There are no imaginary 
parts at > 0. 
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kiQ + UJ k2o + uj 
Figure 20: 2-loop HQET propagator diagram / 




rii n2 



Figure 21: Trivial case 71.5 = 




Figure 22: Trivial case ni = 



^4 7X5 77-4 




Figure 23: Trivial case 713 = 
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If 77.5 = 0, our diagram is just the product of two 1-loop diagrams (Fig. 1^ : 

/(ni, n2, ng, 77,4, 0) = /(rii, ?i3)J(?i2, n^) . (3.12) 

If ni = 0, the inner loop gives G(n3, n5)/(-/ci)"3+"5-'^/2, and fFig. 1^ 

/(O, na, 77,3, 77,4, 775) = G(?73, 775)1(772, 7Z4 + ri3 + 7^5 - rf/2) . (3.13) 

If Us = 0, the inner loop gives /(tzi, 775)/(— A;o)"^'^^"^~'^, and (Fig. |221) 

/(77i, 7Z2, 0, 774, TZs) = J(77i, 715)1(712 + 77i + 27Z5 - C?, 774) . (3-14) 

The cases 772 = 0, 714 = are symmetric. 

When all Ui > 0, we shall use integration by parts The differential operator d/dk2 
applied to the integrand of ()3.11|1 acts as 

^--^^ + ^2*. + ^2«=.-W (3.15) 

Ok2 L>2 UJ JJ4 1^5 

Applying {d/dk2) ■ /c2, {d/dk2) ■ {k2 — ki) to the integrand of ()3.1H) . we get vanishing 
integrals. On the other hand, using k2V /u = D2 — 1, 2{k2 — ki) ■ k2 = — D4 — D^, we 
see that these differential operators are equivalent to inserting 

d - 7Z2 - 7Z5 - 27Z4 + ^ + ^{D^ - D4) , 
d - 772 - ?74 - 2?75 + J^Di + 7^(^*3 " -D5) 

JJ2 

under the integral sign. We obtain the recurrence relations 

[d-n2-n^- 2^4 + 7232+ + n55+(3~ - 4")] / = , (3.16) 

[rf - 712 - ^4 - 2715 + ^22+1" + 7Z44+(3" - 5")] / = (3.17) 

(two more relations are obtained by (1 <-> 3, 2 4)). Similarly, applying the differential 
operator {d/dk2) ■ v is equivalent to inserting 

-2^ + ^Auj\D2 - 1) + ^Au\D2 - D,) , 

LJ2 J^A -t^5 

and we obtain 

[-27122+ + n44+(2" - 1) + n55+(2" - 1")] / = . (3.18) 

(there is also the symmetric relation, of course). 

We can obtain a relation from homogeneity of the integral ()3.11|) in uj. Applying the 
operator ui ■ d/dui to a;~"i~"'2/ p.nj) . -we see that it is equivalent to the factor 2{d — — 
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77.4 — n^) — Til — n2. On the other hand, exphcit differentiation of {—uDi) ^^{~u}D2) 
gives —rii/Di — n^jD^- Therefore, 



\l{d — — — ns) — rii — n2 + nil^ + ^22^] J = . 



(3.19) 



This is nothing but the sum of the {d/dk2) ■ k2 relation (j3.16p and its mirror-symmetric 
{d/dki) ■ ki relation. A useful relation can be obtained by subtracting 1~ shifted (I3.19|l 
from (jTTTjl : 



\d — rii — n2 — — 2n^ + 1 — {2{d — — — n^) 
+ n44+(3" - 5^)1/ = 0. 



ni-n2 + l)V 



(3.20) 



Expressing J(ni, 77,2, ns, 724, ns) from ()3.17jl 

n22+l- +n44+(3- - 5-) 



/(ni,n2,n3,n4,n5) 



d — n2 — n^^ — 2n^ 



(3.21) 



we see that the sum rii+n^+n^ reduces by 1 (Fig.|Hl). Therefore, applying ()2.17p sufficiently 
many times, we can reduce an arbitrary / integral with integer indices to a combination of 
integrals with ns = (Fig. [HI (EH), r^i = (Fig. 121 (Em), rig = (Fig. ESI (EH)- Of 
course, if max(r22,Ti4) < max(ni,n3), it may be more efficient to use the relation mirror- 
symmetric to ()3.17|1 . The relation ()3.2()|1 also can be used instead of ()3.17|1 . 
The second topology of 2-loop HQET propagator diagrams is (Fig. | 



d'^ki d'^k2 



kio + to k2o + 



7r^(-2a;)2('^-*-5) J(ni, n2, rig, n^, n,) 

{ki + k2)o + UJ 2 



(3.22) 



LU 



= -kf , = -k'2 



This diagram is symmetric with respect to (1 ^ 2, 4 ^ 5). It vanishes if 71.4 < 0, or < 0, 
or two adjacent heavy indices are non-positive. If n-^ = 0, our diagram is just the product 
of two 1-loop diagrams (Fig. I^IJ. If ni = 0, it is also trivial (Fig. 




Figure 24: 2-loop HQET propagator diagram J 
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The denominators in ()3.22j) are linearly dependent^: Di-\- D2 — D^, = 1. Therefore pi] 



J 



3~)J. 



(3.23) 



Applying this relation sufficiently many times, we can kill one of the lines 1, 2, 3, and thus 
reduce any integral J with integer indices to trivial cases. In fact, we have not enough 
independent denominators to express all scalar products in the numerator. Therefore, we 
have to consider a more general integral than ()3.22|1 . containing powers of ki ■ /c2 in the 
numerator. This wider class of integrals can also be reduced to the same trivial cases [22] • 
Let's summarize. There are 2 generic topologies of 2-loop HQET propagator diagrams 
(Fig. 123). The integrals with all integer rij can be expressed as linear combinations of 2 
basis integrals (Fig. [221); coefficients are rational functions of d. Here are ra-loop HQET 
sunset diagrams (Fig. I7f[l : 



n{d-2)) 



(3.24) 



2n 





Figure 25: Topologies of 2-loop HQET propagator diagram 




Figure 26: Basis diagrams (all indices equal 1) 




Figure 27: n-loop HQET sunset diagrams 



^HQET denominators are linear in momenta; there are only 2 loop momenta, and 3 denominators 
cannot be independent. 
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3.4 3 loops 

There are 10 generic topologies of 3- loop HQET propagator diagrams (Fig. l28|l . Diagrams 
in the first 2 rows of the figure have one scalar product which cannot be expressed via 
denominators; those in the 3-rd row have 1 linear relation among heavy denominators, and 
hence 2 independent scalar products in the numerator; those in the last row have 2 relations 
among heavy denominators, and hence 3 independent scalar products in the numerator. 
Integration- by-parts recurrence relations for these diagrams have been investigated in [22] • 
They can be used to reduce all integrals of Fig. |2H1 with arbitrary integer powers of de- 
nominators and arbitrary numerators, to linear combinations of 8 basis integrals (Fig. 1201) • 
This algorithm has been implemented in the REDUCE package Grinder j22]^- Recently, 
the heavy-quark propagator and the heavy-light current anomalous dimension have been 
calculated j2S] at 3 loops, using this package. 




Figure 28: Topologies of 3- loop HQET propagator diagrams 

The first 5 basis integrals are trivial: they are expressed via J„ ()3.24|1 (Fig. l27|l and 
Gn fl2.18|l (Fig. ITT|l . and hence via F functions. The integrals in the 3-rd row differ from 
the one in the 2-nd row by replacing a two-loop subdiagram. In the first case, the 2-loop 

^The hep-ph version of this paper contains some corrections as compared to the journal version. 
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I! 






1?.— 





1,1,1,5) 



/lJ(l,l,-l + 2£,l,l) 



Figure 29: Basis diagrams (all indices equal 1, no numerators) 



HQET subdiagram (Fig. |2SI) is substituted: h/i- 



If/{-ujY~'^'^. In the second 



case, the 2-loop massless subdiagram (Fig. ini) is substituted: G2/{—k'^)^ G\/ {—kl^Y 
Therefore, these diagrams are obtained from J3 by multiplying by* 

/2/(6-2c/,l) 3d -7 If 



/2/(5 - 2d, 1) 
^2/(1, 3 -rf) 



2d - 5 /2 ' 
~ d-'i~Gi 



Two basis integrals are proportional to /(1, 1, 1, 1, n) (Sect.ES)) and J(l, 1, n, 1, 1) fSect. lT3j) . 
for non-integer n. The last one is truly 3-loop and most difficult; it will be discussed in 
Sect. 01 



3.5 J(l,l,n, 1,1) 

This diagram (Fig. I30|l has been calculated in [22]. In x-space, it is is 



dt^ dt2 {l{t2 - h)Y"'{lt2Y~\l{t - h)) 



\2-d 



0<tL<t2<t 



'Grinder uses B4 = /3/1//2 and Br, = Ij,G\/G2 as elements of its basis, not the diagrams in Fig. 1291 
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or, going to Euclidean space (tj = —itEi), 

J dh dt2 (t2 - h)''-Hl-'^{t - hf-'^ = Jt— 2^+5 



0<ii<t2<i 



where the power of t is obvious from counting dimensions. Collecting factors from Fourier 
transforms, we have 



T{n-2d + Q)T\dl2-l) ^ 

J(l, 1, 1) = — J , 

Tin) 

where the dimensionless integral J is 

J= j dtidt2{t2-tiY-Hl~\i-hf^ 

0<tl<t2<l 

The substitution ti = xt2 gives 
1 1 





(3.25) 




Figure 30: J(l, 1, n, 1, 1) in x-space 

Integrations in t and x would factor, if not the factor (1 — xt)'^~'^. Let's expand this 
troublesome factor according to the Newton binomial: 



'i-xtf-'=Y: 



k=0 



Here the Pochhammer symbol is 



fc-i 



l[{x + 



i=0 



T{x + k) 
Tix) 



(3.26) 



Then the integrations become trivial: 

{d-2)k 



J = V{n)Y, 



{d-2)k{n-d + 3)k 

,n-d + k + 3)r(n + k + 1) n{n - d + 3) {n + l)k{n -d + A)^ 
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(3.27) 



Recalling the definition of the hypergeometric function 



3-^2 



ai, 02, 03 



X 



fe=0 



A:! 



(6i)fe(&: 



•2 k 



(3.28) 



and taking into account (1)^ = fc!, we can rewrite the result as 



J 



n{n — (i + 3) 
and hence ()3.25j) 



3-^2 



— 2,n — (i + 3 
n+l,n — ci + 4 



Ji^^^^^^^^^)- (n-rf + 3)r(n + l) n+l,n-rf + 4 



(3.29) 



This is the easiest diagram calculation in the world having 3F2 as its result: just a double 
integral. 

It is not more difficult to obtain the general result^ 



J(ni,n2,n3,n4,n5) = 

r(ni + n2 + Us + 2(77.4 + — d))T{ni + ns + 2n4 — d)T{d/2 — ni)T[d/2 — n^) 
r(ni +n2 + n^ + 2n4 - d)T{ni)T{n^)T{ni + 723) 

Til, d — 2n5, rzi + n3 + 2?74 — d 
rii + n3, rii + 772 + n3 + 2n4 - d 



(3.30) 



X 3^2 



1 



In order to calculate the 2-loop diagram of Fig. [23 with a 1-loop insertion in the middle 
heavy line, we need J(l, 1, + 2^, 1, 1). It is easy to shift 77.3 by ±1 using ()3.23j) . Therefore, 
it is sufficient to find it just for one n. The simplest choice for which the algorithm of 
e-expansion (Sect. ISJ works is = 2. We obtain [22] from ()3.29p 

J(l,l,2 + 2£,1,1) ^ 



X 



3{d-4:){d-5){d-6){2d-9) 
e] „ / l,2-2e, l + 4£ 
3 



T{l + 2e) ^ \ 3 + 2£,2 
Expansion of this 3F2 in e is (Sect. ^ 



1 



(3.3i; 



3-^2 



(3.32) 



l,2-2£,l + 4£ \ 
3 + 25,2 + 45 1 J = 2 + 6(-2C2 + 3)5 

+ 12(10C3 - IIC2 + 6)52 + 24(-28C4 + 55C3 - 27(2 + 9)5^ 
+ 48(94(5 - I6C2C3 - 154C4 + 135C3 - 45(2 + 12)5^ + ■ ■ • 

It is not difficult to find several additional terms. 

^I've presented only the result for ni = 112 = 1 in I've no idea why: the general case is not more 
difficult. 
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3.6 /(l,l,l,l,n) 

This diagram (Fig. l3T| has been calculated in j23] , using Gegenbauer polynomials in x- 
space [Hj. This is, probably, the simplest example of applying this useful technique, 
therefore, we shall consider some details of this calculation. In Euclidean x-space, this 
diagram is 

* r 1 97r'^/2 
dU / d'^X = = — - Jt2(n-d)+5 

where \ = d/2 — 1, and the power of t is obvious from counting dimensions. The middle 
vertex has been chosen as the origin, because the propagator with the non-standard power 
begins at this vertex (it is much easier to handle propagators beginning at the origin). 
Collecting factors from Fourier transforms, we have 

jl, , , , ^ 2 T{2{n-d + ?,))T{d/2-n)T\d/2-l) ^ 
The dimensionless integral / is 

1 oo 

^2n-l 

/ = I dt I dxdx 







[{x - vt)^]^ [{x + v{i ~ t)y]^ ' 



where 

d cc — , ~ ~, T" X dec dec ^ 
T{d/2) 

and the angular integration measure dx is normalized: 
dx = 1 . 

x 




v{ti -t) vti 
Figure 31: /(1, 1, 1, 1, n) in x-space 

This problem is more difficult than the one in Sect. 13.51 now we have d- dimensional 
integration over the coordinates of the light vertex. It can be simplified by expanding 2 
massless propagators (whose ends are not at the origin) in series in Gegenbauer polynomi- 
als [Hj 

oo 

J2T'{x,y)C',{x-y), (3.34) 



[{x-yy]^ max2A(x,?/)^ 
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where 



T{x, y) = min 



X y 
y ' X 



Then the angular integration can be easily done using the orthogonality relation 

dxC^^ia ■ x)Cl{h ■ x) = 5,,k,j^C^,ia ■ h) . (3.35) 

In our particular case, the unit vectors a and h are v and — f , and the result involves 
Now we have a single sum: 



d-2 

T{d-2) 



k=0 



-l)T{d + k-2) 
~k\ d + 2k-2 



hit) 



hit) 



x^"-^ [Tix,t)Tix,l-t)r 
[max(a;, t) max(a;, 1 — t)]"'"^ 



The contributions of the regions t < ^ and t > are equal, and we may consider the first 
of them and double it. The radial integral (in x) has to be calculated in several intervals 
separately, because of max and T in ()3.34|1 : 



„2n- 



hit) = I dx * r, + 

J [t(l-t)]'"' 
t 



dx- 



X 



2n- 



-1 rt X -]> 
Ix i-t\ 



[xil-t)] 



d-2 



+ 



dx- 



X 



2n- 



-1 rt i-ti' 

Lx X \ 



\X^ 



id-2 



The result is 



d-2 



1/2 



(d - 2n - 2)r(rf - 2) 



dt^^^Tid + k-2) 







fc=0 



k\ 



X 



-d+2n+k+2l 



t) 



-d~k+2 



t'^il-t) 



~2d+2n-k+4' 



n + k 



d-n+k-2 



Unfortunately, the integral in t cannot be immediately evaluated. The authors of j21] 
had to do some more juggling to obtain 



/(l,l,l,l,n) 



n — 1 X 



Vi2n - 2rf + 6) / 1, - 2, 2n - + 3 

^ ' n + l,2n-(i + 4 



(2n-(i + 3)r(n 



1 - 



(3.36) 

r(rf-n-2)r2(n-d + 3) 
Tid-2) 
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This result can also be rewritten [21] as 



/(l,l,l,l,n) 



r(f-i)r(f-n-i) 

T{d-2) 



3-^2 



n — d + 3,n — d + 3,2n — 2d + 6 
n- d + A,3n-2d + Q 



r{2n - d + 3)r{2n - 2d + 6) 
{n - d + 3)T{3n - 2d + Q) ^ 

- T{d - n ~ 2)T\n - d + 3) 



(3.37) 



this series has larger region of convergence. 

In order to calculate the 2-loop diagram of Fig. 120 with a 1-loop insertion in the middle 
light line, we need /(1, 1, 1, 1, n + e). It is easy to shift by ±1 using the relation [22] 



[{d - 2n5 - 4)5+ - 2{d -n,- 3)] /(1, 1, 1, 1, ng) = 
[{2d - 2^5 - 7)1~5+ - 3~4+5~ + 1"3+] /(1, 1, 1, 1, rig) 



(3.38) 



which follows from integration-by-parts relations (Sect. 13. 3|) (all terms in its right-hand 
side are trivial). Therefore, it is sufficient to find it just for one n. The simplest choice for 
which the algorithm of e-expansion (Sect.E]) works is to use n = 1 in ()3.37j) . We obtain j2Hj 



1(1,1,1,1,1+5) 



4r(l-e) 
9{d-3){d-Ay 



r(i + 4£)r(i + 6£) 
r(i + 7€) 



3-^2 



35, 3e, 6e 
1 + 35, 1 + 75 



1 -1^(1 + 35)1(1-35) 



(3.39) 



This 3F2 function has 3 upper indices ~ 5; therefore, its expansion (apart from the leading 
1) starts from 0{e^): 



3-^2 



35, 35, 65 

1 + 35, 1 + 75 



1=1 + 54(36^ - 513(^6^ + 54(25(5 + 28(2(3)^^ + ■ ■ ■ (3.40) 



It is not difficult to find several additional terms. 



4 Massive on-shell propagator diagrams 

In this Section, we shall consider self-energy diagrams of a massive particle with mass m, 
when the external momentum is on the mass shell: p"^ = m?, ot p = mv. Why this special 
case is interesting? First of all, in order to calculate any S'-matrix elements, one has to 
consider mass fl^ 122] and wave-function [23 I2H] renormalization in the on-shell scheme. 
Calculations of form factors and their derivatives at the point q = belong to this class, 
e.g., the anomalous magnetic moment [22] and electron charge radius [SU] in QED. Finally, 
such calculations are used for QCD/HQET matching |^ IH^. 
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4.1 1 loop 

The 1-loop on-shell propagator diagram (Fig. l32|l is 
It vanishes if rii < 0. 




k + mv 



Figure 32: 1-loop on-shell propagator diagram 

The definition ()4.ip of M(ni, n2) expressed via the dimensionless Euclidean momentum 
K = kE/m becomes 

"^'^ 7r'/'M{n,,n,). 



The definition ()3.7|1 of the 1-loop HQET integral I{ni,n2) expressed via the dimensionless 
Euclidean momentum K = kE/{—'2uj) is 



;i - 2iKqY^{k^y^ 

Inversion ()2.(ij) transforms the on-shell denominator to the HQET one: 

- ^^K, = (4.2) 

(and vice versa). Therefore, the problem of calculating M(ni, 77-2) reduces to the previously 
solved one for I{ni,n2) (Sect. 



, ^ rf , ^ r - rii - 2^2 r -(i/2 + ni + , . 

M{ni,n2) = I{ni,d-ni-n2) = r (4.3) 

r(ni)r(d - rii - n2) 



(Fig.EHl). 
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d — rii — n2 




Figure 33: Inversion relation 



4.2 2 loops 



There are 2 generic topologies of 2-loop on-shell propagator diagrams. The first one is 
(Fig. El: 

^^^ifh = -vrW('^"S"«)M(ni,n2,n3,n4,n5), (4.4) 



Di = vr? — (ki + mvY , D2 = — {k2 + mvY , 
D^ = -kl, D, = -kl, D^ = -{ki-k2f. 



The second topology is (Fig. 1^ : 

.^rf^2(d-E "0 Ar(ni, ris, rig, ^4, n^) , (4.5) 



i^i = — {ki + mf )^ , D2 = — {k2 + mvY , ^3 = rn? — {ki + k2 + mvY 
= -kl , = -k^ 



2 • 



The integrals of both generic topologies (Fig. IHUjl with all integer can be reduced, 
by using integration-by-parts recurrence relations, to linear combinations of 3 basis inte- 
grals (Fig. inZj); coefficients are rational functions of d. Only the topology involve the 
last, and the most difficult, basis integral in Fig. EZI The reduction algorithm has been 
constructed and implemented as REDUCE package RECURSOR in It also has 

been implemented as FORM package SHELL2 

The first 2 basis integrals are trivial. The n-loop on-shell sunset M„ (Fig. l38|) is 

^ ^ i.nd - 4n + 1)2(^-1) 

" + 1 - 4)n ((^ + l)f - 2" - 1)„ K - 2n + l)n-i {n-{n- l)f 
r(l + in- l)e)r(l + ne)Til - 2ne)r"(l - e) 

T{l-ne)T{l-{n + l)e) ' ^ ' ' 

In the 1-loop case, it can be reduced to the 1-loop massive vacuum diagram 
(Fig. El. 
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ki + mv k2 + mv 



Figure 34: 2-loop on-shell propagator diagram M 




k2 



Figure 35: 2-loop on-shell propagator diagram 




Figure 36: Topologies of 2-loop on-shell propagator diagrams 



Figure 37: Basis diagrams (all indices equal 1) 
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Figure 38: n-loop on-shell sunset diagrams 

ld-2 




2d-3 



Q 



Figure 39: 1-loop on-shell diagram 

Instead of using A^(l, 1, 1, 0, 0) (the last diagram in Fig. EH which is divergent) as an 
element of the basis, one may use the convergent diagram A^(l, 1, 1, 1, 1) (the second one 
in Fig. IH^. It has been calculated in 



iV(l,l,l,l,l) 



4r^(i + £) 

3e{l - 4£2) 
1 



-2 3-F2 



1 i -£ i-e 



2 ^ ' 2 



+ 



1 i i 

^ ' 2 ' 2 



^ 3(1 + 2£) / r(i -4£)r(i + 2£)r^(i _ ^ 



r(i -3£)r(i-2e)r(i + £) 



Its value at £ = is ^Hl 



(4.8) 



iV(l,l,l,l,l) = 7r2log2--C3 + 0(£) 



(4.9) 



D. Broadhurst discovered [22] a symmetry group of a class of 3F2 functions like those 
in ()4.8j) : this allowed him to obtain ()4.9j) and two further terms, up to 0{e'^), by purely 
algebraic methods. 

Using inversion of Euclidean dimensionless momenta ()2.6|) and ()4.2p . we obtain 
relation of Fig. EHl 



^4 






Kd — n2 — 








^5 \ 



n2 



Figure 40: Inversion relation 
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4.3 2 loops, 2 masses 

In realistic theories, there are several massive particles with different masses, as well as 
massless particles (QED with e, r; QCD with c, b, t; the Standard Model). Therefore, 
diagrams like in Fig. [^appear. Here the thin line is massless, the thick solid line has mass 
m, and the thick dashed one has mass m'. Combining identical denominators, we obtain 
the topology shown in Fig. 

Figure 41: 2-loop on-shell propagator diagram 




Figure 42: Topology of 2-loop on-shell propagator diagrams with 2 masses 

An algorithm for calculating such diagrams based on integration-by-parts recurrence 
relations has been constructed in |2Z1- The algorithm has been implemented in REDUCE. 
There is one scalar product which cannot be expressed via the denominators. The in- 
tegrals of Fig. |42l with arbitrary integer powers of the denominators, and with arbitrary 
non-negative power of the scalar product in the numerator, can be expressed as linear 
combinations of 4 basis integrals (Fig. EH|l : coefficients are rational functions of d, m, and 
m'. 




Figure 43: Basis 2-loop on-shell propagator diagrams with 2 masses 

Instead of using the last diagram in Fig.llHlas an element of the basis, one may prefer to 
use the sunset diagram (the third one in I43|l with one of the denominators squared (either 
the one with mass m or the one with mass m'). Therefore, the reduction algorithm for 
the sunset diagrams with the masses m, m', m' and = m?, constructed in [2H1, is just a 
particular case of the algorithm from . 
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The first 2 basis integrals in Fig. are just products of 1-Ioop massive vacuum di- 
agrams, and tlius trivial. Tlie 2 nontrivial ones liave been calculated in [2ZI using tlie 
Mellin-Barnes representation tecfmique: 



/2U-2e 



r2(i + 5) 



+ 



£2(1 -e) \l-2e 

2 \ l-£ 

3-^2 



3-^2 



l,l,-l + 2e 



m 



m' 



12 



m 



h 



m'2 

/2\-2e 



I e - - e 
3 - 2e, I 

1 



m 



m 



12 



r2(i + £) 



2(l-£)(l + 2£) 



3-^2 



3 

-' 2 



m 



2\ i-e 



m' 



12 



l-2e 



3-^2 





1 

2 








2£, 


3 
2 


m'2y 1 



(4.10) 



If m' = m, they are not independent - Ji can be reduced to Jq (the last diagram in 
Fig- EH) using integration by parts (Sect. 



3d -8 



3{d-2y 



(4.11) 



4(rf-4) " 8(d-3)(rf-4) 
We can express A^(l,l,l,l,l) via Jq or Ji fl4.1()|l : in the second case, for example, we have 



iV(l,l,l,l,l) 



+ 



(rf-4)2 
1 



rf-3 



3-^2 



(3rf- 10) 

I e ^-e 
2-2e,f 



1 



_(rf-2)(rf-5) 
1 



3-^2 



1 2f - 

2 - e, I + £ 



(c/-3)(c/-4) 



^ r2(i-£)r(i + 2e)r(i-4e) _ ■ 
r(i + e)r(i-2e)r(i-35) 



(4.12) 



This should be equivalent to 

The basis integrals ()4.10|) . expanded j^ZI in e up to (9(1), 



12(1+5) 



r2(i + 5) 



— m 



2-4e 



^ + 1 + 2(1 -rY(^-^ + L-)-21ogV + ^ 



t2-i£ 



13, 
^ H 21ogr + 6 



0{e) 



m 



-4e 



15 19 
— + - + 2(1 + r)2L+ + 2(1 - r)2L„ - 2 log^ r + y 



+ 0(5) 



(4.13) 
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(where r = m'/m), are expressed via dilogarithms: 



Li2(— r) + i log^ r — logr log(l 



1, 



Li2(-r ^) + logr Mog(l + 
Li2(l-r) + |logV + V 



r 



— Li2(r) + I log^ r — logr log(l — r) + |7r^ (r < 1) 
Li2(r~"'^) + logr""*^ log(l — r""*^) (r > 1) , 



L+ + L_ = i Li2(l - r^) + log' r + j^vr' 



-iLi.fl-r^'N 



12' 



(4.14) 



The sunset basis integral Jq fl4.1()|l (the third one in Fig. EHl) was recently expanded jnOj 
up to e^; the coefficients involve harmonic polylogarithms of r^^. In particular, at r = 1, 
one more term in ()4.9p {0{e^)) can be obtained. 



4.4 3 loops 

There are 11 generic topologies of 3-loop massive on-shell propagator diagrams fFig. I44p. 
Ten of them are the same as in HQET (Fig. OHj) . plus one additional topology with a 
heavy loop (in HQET, diagrams with heavy-quark loops vanish). Integration- by-parts 
recurrence relations for these diagrams have been investigated in |2H] • They can be used to 
reduce all integrals of Fig. |3U with arbitrary integer powers of denominators and arbitrary 
numerators, to linear combinations of the basis integrals in Fig.HSl This algorithm has been 
implemented in the FORM package SHELLS |2HI- Some of the basis integrals are trivial. 
Others were found, to some orders in e, during many years, in the course of calculation 
of the anomalous magnetic moment in QED at 3 loops [SH] • Some additional terms of e 
expansions were obtained in |28j . 

Performing inversion ()2.(ij) of the loop momenta, we obtain the relations in Fig. IIHI 
For example, the HQET ladder diagram with all indices rij = 1 is convergent; its value at 
d = 4 is related [10] to a massive on-shell diagram fFig.l47|) by the second inversion relation. 
This is one of the basis diagrams of Fig. and its value at d = 4 is known (Fig. 
Calculating this ladder diagram with Grinder (Sect. I3.4|l and solving for the most difficult 
HQET 3-loop basis integral in Fig. EHl we obtain the e expansion of this integral up to 
0{e). This concludes the investigation of the basis integrals of Fig. 123 (Sects. I3.4H3.6|) . 
and allows one to solve 3-loop propagator problems in HQET up to terms 0{1) (see [22] )• 
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Figure 44: Topologies of 3-loop massive on-shell propagator diagrams 

<Z> ^ O ^ O <l> 
<^ ^ ^ <^ ^ 

© ^ ^ ^ W <s> 

Figure 45: Basis diagrams (all indices equal 1, no numerators) 



ni n2 




d — — Uq — Tij — rig 




d — rii — V *Kd — n2 — 

- ^6/ "-6 n7 \ _ m 



rii n-s n2 



Figure 46: Inversion relations 




5C5 + I2C2C3 



Figure 47: Inversion relation: all = 1, d = 4 



5 Hypergeometric functions and multiple ( values 

In this Section, we shall consider some mathematical methods useful for expanding massless 
and HQET diagrams in e. 

5.1 Multiple ( values 

As you may have noticed, expansions of various Feynman integrals in e contain lots of 
C-functions. They appear when expanding F-functions: 



They also appear, together with their generalizations - multiple ( values, when expand- 
ing hypergeometric functions of unit argument with indices going to integers at e — > 
(Sect. 15. 2| I5.3jl . Such hypergeometric functions appear in massless and HQET diagrams 
(Sect. 13.5^ imjl . In this Section, we shall briefly discuss multiple ( values; more details 
can be found in jl^, where more general functions are also considered. 
The Riemann ^-function is defined by 




00 



r(l + e) = exp 



(5.1) 



n=2 




(5.2) 



n>0 
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Let's define multiple ( values as 

c = T ^ 

SsiS2 / , si 



ni>n2>0 



S2 ' 



1 "-2 



(5.3) 



ni>n2>n3>0 



Hi n2 



and so on. These series converge at si > 1. The number of summations is called depth k; 
s = si + ■ ■ ■ + Sfc is called weight. All identities we shall consider relate terms of the same 
weight (where the weight of a product is the sum of the factors' weights). 
For example, there is 1 convergent sum at weight 2: 

C2, 

2 at weight 3: 

Ca , C21 , 
4 at weight 4: 

C4 5 Csi 7 C22 5 C211 , 
and 8 at weight 5: 

Cs , C41 , C32 , C23 , Csii , C221 , C212 , C2111 • 
Suppose we want to multiply (s(siS2- 

ni>n2>0 

Here n can be anywhere with respect to rii, 77-2. There are 5 contributions: 

/ . Q Si S2 S,SS\S2 1 



n>ni>Ti2>0 



V ' -c 



n°7i^ rir, 

n=ni>n2>0 



S2 5 



1 "-2 



S2 SS1SS2 ) 



y — 

ni>n>n2>0 

/ J e St So SSI, 

n=n2>0 ^ ^ 



ni>n=n2>0 



S2 ) 



ni>n2>n>0 



Wn^ n. 



1 '"2 
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Therefore, 



CsCsiS2 CssiS2 ~l~ Cs+si,S2 Csiss2 ~l~ Csi,s+S2 ~^ Csi 



S2S 



(5.4) 



This process reminds shuffling cards fFig. I48|). The order of cards in the upper deck, 
as well as in the lower one, is kept fixed. We sum over all possible shufflings. Unlike real 
playing cards, however, two cards may be exactly on top of each other. In this case, they 
are stuffed together: a single card (which is their sum) appears in the resulting deck. A 
mathematical jargon term for such shuffling with (possible) stufflng is stujfling, see |41j . 




Figure 48: Stuffling: shuffling with (possible) stufflng 
For example, the stuffling relations for weights s < 5 are 



(^2 



C22 + C4 + C22 , 

C2C3 = C23 + Cs + C32 , 

C2C2I = C22I + C4I + C22I + C23 



(5.5) 



,212 
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Now we are going to derive an integral representation of multiple ( values. Let's consider 
the integral 



1 Xl X2 X3 

dxi f dx2 f dxs f dx4 



Xl J X2 J X3 J X4 


1 XI X2 

dxi f dx2 f dxs ^ 



■x^ 



Xl J X2 J xs n 
000 

1 Xl 

dxi f dx2 „, 1 _ 



■^2 " 9 
Xl J X2 ^ 



1 

Xl ■ — 

Xl n-^ 



1 

^' 

It is easy to guess that 

1 f dxi dxs 



-Xs 

n" J Xl Xg 

l>Xl>--->Xs>0 



, dxi dxs-i dxs 
Cs= / • (5.6) 



But we need the sum ()5.2j) . This is also easy: 

dxi dxs-i dxs 

Xl Xs-l I- X, 

l>a;i>->2,'s>0 

Let's introduce short notations: 
dx dx 



LUo= , UJi = • 

X 1 — X 

All integrals will always have the integration region 1 > > • • ■ > a;^ > 0. Then (j5.6|l 
becomes 



Cs = j o^r'^i • (5.7) 
This can be generalized to multiple sums: 



CsiS2S3 — I ^^1 ^^1 '^^1 , 



(5.8) 
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and so on. 

Suppose we want to multiply (2(2- 



C2 = J ^0^1 • J ^0^1 ■ 

l>a;i>X2>0 l>x[>x'2>0 

The ordering of primed and non-primed integration variables is not fixed. There are 6 
contributions: 



■,22 5 
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,31 ; 



1 > Xi > X2 > x'^ > X2 > : J UqUiUquJi = (2 
1 > Xi > > X2 > X2 > : J cjo^^o^^it^i = C31 

1 > Xi > x']^ > X2 > X2 > : J LdQUJoUJiUi 

1 > x'^ > Xi > X2 > X2 > : J cuqujqujiUJi = C3 

1 > x'l > Xi > X2 > X2 > : J cuqujoujiUi = C31 , 

1 > x']^ > X2 > Xi > X2 > : J laJqLaJicuoUi = C22 ■ 
Therefore, 

a = 4C31 + 2C22 . (5.9) 

Now we are multiplying integrals, not sums. Therefore, our cards are now infinitely 
thin, and cannot be exactly on top of each other. There are just 2 kinds of cards: uq and 
Ui, and we sum over all possible shufflings of two decks (Fig. 

For example, the shuffling relations for weight 5 are 

C2C3 = 6C3II + 3C22I + C2I2 , 

C2C2I = C4I + 2C32 + 2C23 ■ 

The integral representation allows us to derive another set of useful relations, even 
simpler than shuffling - duality relations. Let's make the substitution Xj — 1 — Xj. Then 
ujq ^ uji, in order to preserve the ordering 1 > Xi > ■ ■ ■ x^ > 0, we have to arrange 
all the 00 factors in the opposite order. In other words, after writing down an integral 
representation for a multiple ( value, we may read it in the Arabic fashion, from right to 
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Figure 49: Shuffling 
left, simultaneously replacing uq ^ ui. 
Cs = y" ^^o^o^i = J ujquJiUJi = C21 , 
C4 = y t^o^^ot^ot^i = j ujqUiujiuji = C211 , 
Cs = y uJoUJoUJoUJoUl = J u;ot^it^it^it^i = C2111 , 
C41 = j ujqujoUqUiUi = j u;ou;ou;iu;iu;i = Csii , 
C32 = J ujqujoujiujo^i = J ujquJiUJoUJiUJi = C221 , 
C23 = J ujqujiujo^^o^^i = J ujQUJiUJiUJoUJi = C212 • 
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Duality relations are the only known relations which say that two distinct multiple ( values 
are just equal to each other. 

Let's summarize. There are 3 distinct multiple ( values of weight 4, 

C4 = C211 5 C31 ) C22 5 
due to duality ()5.11|) . The stuffling ()5.5|) and shuffling ()5.9p relations yield 

Cl = C4 + 2C22 C22 

CI = 4(31 + 2(22 ^ Cai 
where we have used 

This follows from 

^'-y ^'-90- 

Therefore, all multiple ( values of weight 4 can be expressed via (4. 
There are 4 distinct multiple ( values of weight 5, 

Cs = C2111 5 C41 = Can ) C32 = C221 5 C23 = C212 5 
due to duality ()5.11|) . The stuffling ()5.5|) and shuffling ()5.10|) relations yield 

C2C3 = C32 + C23 + CS ! 

C2C3 = 6C41 + 3^32 + C23 , 
C2C3 = C2C21 = C41 + 2C32 + 2C23 • 

Solving this linear system, we obtain 

C41 = 2C5-C2C3, C32 = -yC5 + 3C2C3, (23 = ^(5-2(2(3- (5.13) 

Therefore, all multiple ( values of weight 5 can be expressed via (5 and (2(3- 

Weight 5 is sufflcient for calculating finite parts of 3-loop propagator diagrams. We 
shall not discuss higher weights here. 

5.2 Expanding hypergeometric functions in e: example 

There is an algorithm to expand in e hypergeometric functions of unit argument whose 
indices tend to integers at e — >■ 0. The coefficients are linear combinations of multiple 
( values. We shall first consider an example: expanding the 3F2 function (j2.2H|l . which 
appears in 1, 1, 1, 2 + e) (Sect. 1231), up to 0{e). 



3. 



(5.12) 
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The function we want to expand is 

(2-2£)„(2 + 2£)„ 
(3 + 5)„(3 + 2£)„ • 



n=0 



When arguments of two Pochhammer symbols differ by an integer, their ratio can be 
simphfied: 

(2 + 2e)n 2 + 2e 



(3 + 2e)n n + 2 + 2e' 
Step 1. We rewrite all Pochhammer symbols (m + le)n via (1 + le)n+m-i- 

(2 - 2£)„ = ^if^^ , (3 + e) - + 



Then our function F ()5.14|) becomes 

l-2£ ^ -r -r ; . Z^„ + 2 + 2£ l+£„+2 ' ' 

n=U 

In what follows, we shall consider F', and return to F at the end of calculation. It is 
convenient to introduce the function Pn{s)'- 

n 

(1 + e)n = n\ P„(£) , P„(£) = J] (^1 + i-) . (5.16) 

n'=l 

Then 



1 Pn+l(-2g) 

n + 2)(r2 + 2 + 2e) P„+2(e) 



n=0 

Step 2. We expand the rational function in front of P's in e: 
1 1 2£ 

+ ■■■ 



(n + 2)(n + 2 + 2e) (n + 2)2 (n + 2)^ 

In this simple case, all n-dependent brackets are the same (n + 2); in more general cases, 
we should have to decompose each term into partial fractions with respect to n. Then 

P' - 1 Pn.+l(-2£) _ ^ 1 , a^/^2n 

^ - + 2)2 P„^,(.) ^ + 2)3 + ^ ' 

because we may replace all Pn{l£) ^ 1 in the 0{e) term. Shifting the summation indices, 
we get 



P„(£:) ^ n3 

n=2 ^ ' n=2 
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Step 3. We rewrite all Pn+mi^^) via P„_i(te): 

P„(£)= (1 + ^)P„_1(£). 

We again expand in e (and, if necessary, decompose into partial fractions) rational functions 
in front of P's: 

n=2 ^ ' n=2 

Step 4. We add and subtract a few first terms to have all summations start from 
n = 1: 

^'-|:i^t^-i-3(C3-i)..o(.Y 

Step 5. We expand P„_i(fe) in e: 

Pn^x{e)= J] (l + l^) = l + ^i(n)£ + ^n(n)£' + --- (5.17) 

n>n'>0 

where 



S2 ' 



n>ni>n2>0 ^ ^ 

and so on. These finite z sums obey the same stuffling relations as infinite C sums, for 
example, 

Zs{n)z^^^^{n) = Zssjs^in) + Zs+si,s2i'n') + Zs.ssiin) + Zs^^s+S2{n) + Zs,s2s{n) ■ (5.19) 

Therefore, coefficients of expansion of ratios of products of P functions can always be made 
linear in z sums. In our simple case, 

Pn-i{-2e) l-2z^{n)e + - 



Pn-i{e) l + zi{n)e 
Therefore, 



l-3zAn)e + 



l)e + - 



n>0 n>ni>0 



= C2-1-3(C21 + C3-1)^+--- 

Recalling the duality relation ()5.1ip (21 = Cs? "we obtain for our function F ()5.14|1 

P = 2(2 + 96 + ■ ■ ■ ) [C2 - 1 - 3(2C3 - 1)£ + ■ ■ ■ ] 
= 4(C2-l) + 6(-4C3 + 3C2-l)£ + --- 
Thus we have reproduced the first two terms of ()2.23p . 
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5.3 Expanding hyper geometric functions in e: algorithm 



Now we shall formulate the algorithm in a general setting. We want to expand 



oo 



]\i{mi + lie)n 



(5.20) 



n=0 



n 



m e. 



Step 1. We rewrite all Pochhammer symbols (m + le)n via (1 + le)n+ra-i'- 



(m + le) 




n 



(1 + /e) ■ ■ ■ (m - 1 + ■ 



Introducing the function P we rewrite F fl5.2()|l as 



oo 



Hi Pn+mi-lik^) 



Hi' Pn+m'-l{l'ii£) ' 



n=0 



where i? is a rational function of e and n. 
Step 2. We expand 

i?(n, e) = i?oH + Ri{n)e H 

and decompose each Rj{n) into partial fractions with respect to n. There is a catch here: 
partial-fractioning of a convergent series can split it into a combination of logarithmically 
divergent ones. Therefore, it is necessary to keep the upper limit of summation finite, 
and go to the limit only after combining such series together. Because of shifts of the 
summation variable, cancellation of divergent series leaves us with a few terms around the 
upper limit. They, however, tend to zero, and this means that we may formally manipulate 
logarithmically divergent series as if they were convergent. 

Shifting summation indices, we can write F ()5.2()j) as a sum of terms of the form 




Step 3. We rewrite all P„+m(/£:) via P„_i(fe): 




Our F becomes a sum of terms of the form 




oo 
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where i? is a rational function of e and n. Moving from terms with lower powers of e to 
ones with higher powers, we expand i?'s in e: 



R{n,e) = 1 + Ri{n)e H 

Terms of higher orders in e are decomposed in partial fractions again. Finally, F becomes 
a sum of terms of the form 



1 U.Pn-iihe) 



n=no 



Step 4. Adding and subtracting terms with n from 1 to no — 1, we rewrite F as a sum 
of terms of the form 

and rational functions of e (which can be trivially expanded in e). 

Step 5. We expand each Pn-i{ls) as ()5.17|) . Expansions of ratios of products of P 
functions contain products of z sums; they are reduced to linear combinations of z sums 
by stufHing relations (e.g., ()5.19|l ). Using 

oo ^ 
n=l 

we can express coefficients of e expansion of F as linear combinations of multiple C, values 
(discussed in Sect. 15.111 and rational numbers. Q. E. D. 

A few historical comments. In 2000, I needed to expand some 3F2 functions of unit 
argument in e (in fact, they were ()3.32|) and ()3.40|) ). 1 asked David Broadhurst how 
to do this, and he replied: just expand Pochhammer symbols in e, the coefficients will 
be expressible via multiple C, values. Following this advise, I implemented the algorithm 
described above in REDUCE, and obtained the expansions ()3.32|1 and ()3.40|1 . Presentation 
in Sects. ini21 and EiSl closely follows my notes from 2000. Later this algorithm was described 
as Algorithm A in jl2I (this paper also contains more difficult algorithms (B, C, D) to 
expand some more complicated sums in e). All of these algorithms are implemented in the 
C++ library nestedsums jl^l^^- It is based on the C++ computer algebra library GiNaC |44j . 
In order to use nestedsums with recent versions of GiNaC, you should download a recent 
version from the URL shown in 

Unfortunately, no general algorithm is known to expand in e hypergeometric functions 
of unit argument some of whose indices tend to half-integers at e — 0. Such functions 
appear in on-shell propagator calculations, see Sects. 14.21 14.31 Known expansions show 
that many constants more complicated than multiple C, values appear. 

^^Unfortunately, there exist different notations for multiple C, values. We follow here; in 02 ESI the 
order of indices is reversed. 
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